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We present full atomistic calculations of the spin- flip time (Ti) of electrons and holes mediated by 
acoustic phonons in self-assembled Ini_a;Gaa;As/GaAs quantum dots at zero magnetic field. At low 
magnetic field, the first-order process is suppressed, and the second-order process becomes dominant. 
We find that the spin-phonon-interaction induced spin relaxation time is 40 - 80 s for electrons, and 
1-20 ms for holes at 4.2 K. The calculated hole-spin relaxation times are in good agreement with 
recent experiments, which suggests that the two-phonon process is the main relaxation mechanism 
for hole-spin relaxation in the self-assembled quantum dots at zero field. We further clarify the 
structural and alloy composition effects on the spin relaxation in the quantum dots. 
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The electron/hole spin in semiconductor quantum dots 
(QDs) is believed to be a promising candidate for solid 
state quantum computations [l[. Recently years, huge 
progress have been made experimentally in initiahzation, 
manipulation and controUing of the spins in QDs0-Q. 
However, the short spin hfetime is stih a major obsta- 
cle to realize the quantum computation. Until now the 
mechanisms for spin relaxation in QDs are still not well 
understood both in theory and in experiments. Spin- 
phonon interaction due to the spin-orbit coupling (SOC) 
effects is one of the main mechanisms lead to spin relax- 
ation in the QDs d, 01 , which depends strongly on the 
geometries and compositions of the QDs. An accurate 
description of SOC is crucial to understand the spin re- 
laxation. Unfortunately, the understanding of the spin 
relaxation from the atomistic level is still unavailable. 

Electrons were expected to have very long spin lifetime, 
because they have small SOC in QDs. However, because 
of hyperfine interactions with nuclear spins, the electron 
spin coherence time is greatly reduced ('~ 500 ps) 
The electron spin lifetime can be prolonged by applying 
an external magnetic field to polarize the nuclear spin 
0, Q . Unfortunately, at the same time, the electron spin 
lifetimes (Ti) decrease fast with the magnetic field, being 
proportional to In contrast to electrons, holes 

have p-like atomic wave functions. Therefore, the hyper- 
fine interactions with the nuclear spins are small [8| . In- 
deed, recent experiments show that holes have very long 
spin lifetimes in InAs/GaAs QDs. For example, Heiss et 
al. found that the holes can have lifetimes of 0.27 ms 
at 1.5 T Q. Only very recently, high fidelity spin states 
have been prepared at low magnetic field, taking the ad- 
vantages of the extreme long hole spin lifetimes (■^ 1 ms) 
i- 



At very low magnetic field, the first-order spin-phonon 
interaction is greatly suppressed, and the multi-phonon 
process becomes dominate. The Ti cx law breaks 
down at low magnetic field |3, @|- Most previous theo- 
retical works focus on spin relaxation in a large magnetic 
field and there are only few studies of spin relax- 

ation in QDs at low magnetic field ^ . These studies are 
all based on effective mass approximations or k • p theory 
[S [13, III) in which the effective SOC are added in by 
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hand in the form of Dresselhaus interactions and Rashba 
interactions. These continuum theories treat poorly the 
local strain and alloy composition effects especially for 
hole, which may play an extremely important role to the 
effective SOC. 

In this letter, we present the first atomistic calcula- 
tions of the spin relaxation of electrons and holes at 
zero magnetic field in self-assembled Ini_a;Gaa;As/GaAs 
QDs. In this method, the Hamiltonian naturally takes 
the full SOC effects into account at the atomistic level. 
The method needs no fitting parameters to the QD exper- 
iments, and therefore can give quantitative predictions 
that can be directly compared to experiments. Our cal- 
culations show that the holes have extremely long spin 
lifetimes, 1-20 ms at r=4.2 K, which is in excellent agree- 
ment with the recent experimental data This con- 
firms that the two-phonon scattering process with the 
SOC is the main mechanism for the spin relaxation of 
holes at zero magnetic field. 

At low magnetic field, the first-order process is sup- 
pressed, because the Zeeman splitting between spin-up 
and spin-down states is very small. Therefore, the den- 
sity of phonon states that satisfy energy conservation is 
extremely low. In this case, the spin relaxation via the 
two-phonon processes become dominant: a spin at the 
initial (labeled as i) state, e^, absorbs a phonon of mo- 
mentum q and jumps to an intermediate (s) state, ts- It 
emits a phonon with momentum k and relaxes to the fi- 
nal (/) state, £/, which has an opposite spin of the initial 
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state. The electron and hole spin- flip rate (t~^) from the 
initial state to the final state is given by the second-order 
Fermi's Golden Rule, 
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where only long-wave acoustic phonons are involved in 
the process, Wq = Ci/|q|, and are the sound speeds for 
the v= LA and TA modes. Nq = l/[exp(fiwq/fc_BT)-l] is 
the Bose-Einstein distribution function. The summation 
is over the intermediate states includes all the states ex- 
cept for the initial and final states. The matrix elements 
M^^f are given by: 
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where and \ipf) are the wave functions of initial 
and final state, respectively. Q!(q) is the electron-phonon 
coupling strength. We have considered three electron- 
phonon interaction mechanisms in the QDs [lo| : (i) 
electron-acoustic-phonon interaction due to the deforma- 
tion potential — LADP), (ii) electron-acoustic-phonon 
interaction due to the piezoelectric field for the longi- 
tudinal modes {v = LAPZ), and (iii) electron-acoustic- 
phonon interaction due to the piezoelectric field for the 
transverse modes {v = TAPZ). The parameters are taken 
from Ref. The overall spin relaxation time Ti — 

To calculate Ti, it is crucial to obtain high-qu ality 
single-particle wave functions and energy levels [lOj. 
In this work, we use an atomistic pseudopotential 
method, which has been proven to be very accurate for 
InAs/GaAs QDs [l3-[3]- We consider realistic lens- 
shaped Ini_a:Gaa;As/GaAs QDs embedded in a GaAs 
matrix containing 60x60x60 8-atom unit cells. The dots 
are assumed to grow along the [001] direction, on the top 
of the one- monolayer wetting layers. 

We first relax the dot+box system by minimizing the 
strain energy as a function of the coordinate {R„.q} of 
the a-th atom at the site n for all atoms, using the va- 
lence force field (VFF) method 15|, |l6| . Once we have 
all the atom positions, we obtain the energy levels and 
wavefunctions by solving the single-particle Schrodinger 
equation fvi\ . 
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Va (r) is the screened atomic pseudopotential for the a-th 
element, including a local potential and a non-local SOC 
term |12l |. 



^so,c. = ^IO'5Vz^a(r)L.S(Z|, 
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TABLE L Spin relaxation times of the electrons (Tf ) and the 
holes (Ti) in self-assembled Ini_2,Gaa;As/GaAs QDs at 4.2 
K. All lengths are in nm. is the weight of the minority 
spins of the initial states. 
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where / is the angular momentum (only Z=l is used). L 
is the angular momentum operator, and S is the spin 
operator. The atomic pseudopotentials are fitted to the 
bulk properties. Once the parameters are determined, 
there are no free parameters for modeling the QDs. The 
method naturally includes the Rashba and Dresselhaus 
SOC in a "first-principles" manner. 

The single-particle Schrodinger equation is solved by 
the linear combination of bulk bands(LCBB) method 
[l8| . We use eight bands for both the electrons and holes, 
and take the inter-band coupling (heavy-hole-light-hole 
coupling and the valence-conduction band coupling) into 
account, which is very important for hole/electron spin 
relaxation. A 6x6x 16 k-mesh converges very well the re- 
sults. Due to the SOC, the wave functions mix opposite 
spin components, i.e., ipi = a| t) + PI i)- 

To calculate the matrix elements, Mq, we first calcu- 
late the single-particle wave functions on a sparse k mesh 
following Ref. [l^. We then interpolate the wave func- 
tions into a more dense 49x49x129 fc-mesh to get accu- 
rate results. We sum over 20 intermediate states (s in 
Eq. [1]) , which converge a spin lifetime < 1 s for electrons 
and < 1 ms for holes. Details on how to calculate Mq 
will be published elsewhere. 

In Table HI we list the electron and hole spin lifetimes 
for some typical Ini^ajGa^As/GaAs QDs at r=4.2 K. 
We also list the weight of the minority spin of the initial 
states in the Table. From the Table, it can be seen 
that the electron spin lifetimes due to the spin-phonon 
interaction are extremely long, 40-80 s, because the spin 
mixture is quite small (See Table HI. We find that the 
inter-band coupling has a large contribution to the spin 
relaxation for electrons. Ignoring inter-band coupling 
leads to a spin lifetime that is about three times longer 
for electrons. However, the actual electron-spin lifetimes 
are determined by the hyperfine interaction with nuclear 
spins at low magnetic field. 

The hole-spin lifetimes of typical Ini_a;GaxAs/GaAs 
QDs are 1-20 ms. For alloy dots, the typical hole-spin 
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FIG. 1: (Color online) Spin relaxation times of electrons 
and holes in lens-shaped Ini_a;Gaa; As/GaAs QDs of different 
heights [(a), (c)], and base diameters [(b), (d)]. Height varies 
keeping diameter at 25 nm and the diameter varies keeping 
height at 3.5 nm. Black lines are the results for pure dots, 
whereas the red lines are the results for ally dots with a;=0.5. 



lifetime is a few ms, which is in good agreement with 
recent experimental data ('^ 1 ms) 0]. This suggests that 
two-phonon spin scattering is the main mechanism that 
leads to spin flip for holes in these self-assembled QDs at 
low magnetic field. This finding differs greatly from the 
results given in Ref. 
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where the second-order phonon 
process is not considered. The spin lifetimes of the holes 
are about three orders of magnitudes shorter than those 
of the electrons. This is because the single-particle energy 
spacing of the electrons is 3-5 times larger than that of 
holes, and the spin mixing due to the SOC for holes is 
also much stronger than that of electrons, as shown in 
Table m The spin- flip times in pure QDs are significantly 
longer than those of alloy QDs. The main reason for the 
longer hole spin-flip time in pure QDs is that the energy 
spacing in pure QDs is larger, and the spin mixture in 
alloy dots tends to be slightly larger than that in the pure 
dots. 

The atomistic theory reveals some important features 
that can not be observed hy k-p method, such as geome- 
try and alloy compositions effects to the spin relaxation. 
To explore the effects of geometry on the spin flip time in 
QDs, we calculate the spin flip time of lens-shaped QDs, 
as functions of dot height (diameter) while fixing the dot 
diameter (height). The results are presented in Fig. [T] 
(a) and (b) and for electrons, and Fig. [T] (c) and (d) for 
holes. For pure dots (black dots), the spin lifetime of the 
electrons decreases as dot height increases from 2 nm to 
3.5 nm, and becomes flat for taller dots. The spin life- 
times of holes, also decreases with increases in the height 
of the dots. For extremely tall pure InAs/GaAs quan- 
tum dots, the spin flip time is extremely fast (about two 
orders of magnitudes faster). This is because in very tall 
pure InAs/GaAs QDs, holes tend to localize in the in- 
terface between the dots and matrix, resulting in a very 
small energy spacing, and a fast spin relaxation. The spin 
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FIG. 2: Spin relaxation time of Ino.5Gao.5As/GaAs QDs with 
different Ga distributions. 



lifetimes decrease with increasing dot diameter for elec- 
trons [Fig. [11(c)]. In contrast, the spin lifetimes increase 
with increasing dot diameter for holes [Fig. [T] (d)]. To 
understand the effects of geometry on the spin lifetime, 
we note that there are two competing factors: the energy 
spacing and spin mixing. For electrons, the SOC is not 
very sensitive to the dot geometry, and the spin lifetimes 
are mainly determined by the single-particle energy spac- 
ing. For holes, increasing dot height leads to a smaller 
energy spacing between the levels, and the spin mixing 
also increases with dot height because of the stronger 
effective SOC. Therefore the spin lifetime decreases. In- 
creasing the dot diameter, however, leads to a smaller en- 
ergy spacing between the levels, which tends to increase 
the spin flip time, but the spin mixing decreases with in- 
creasing dot diameter. The overall effect is that the spin 
lifetime increases with increasing dot diameter. There- 
fore, for pure InAs/GaAs QDs, a flatter QD has a much 
longer lifetime than a tall QD. 

We also show that the geometry-dependent spin life- 
time for Ino.5Gao.5As/GaAs alloy QDs (red dots) in Fig. 
[1] For electrons, the spin lifetimes of alloy QDs are very 
similar to those of pure QDs, whereas for holes, the alloy 
QDs show a much weaker geometry dependence for their 
spin lifetime. When the dot height increases, the lifetime 
decreases slightly with increasing dot height. However, 
for most of the dots, the alloy dots have shorter lifetimes 
than those of the pure dots. For very tall dots, the alloy 
dots have longer lifetime, because the holes are not lo- 
calized at the dot-matrix interface in the alloy dots. The 
spin lifetimes of holes are also flat with increasing of the 
dot diameter. 

It has been shown by Mlinar and Zunger [2^ that the 
randomness of In, Ga alloys may have large effects on 
their optical properties for finite sizes Ini_a;Gaa;As/GaAs 
QD (< 10^ atoms). To investigate the alloy randomness 
effects on the spin relaxation time, we calculate Ti for 
Ino.5Gao.5As/GaAs QDs, with various random distribu- 
tions of Ga atoms in the dots. The results are shown in 
Fig. [2j We find that the spin lifetimes of electrons are 
distributed between 44 s to 50 s, in which the lifetime 
fluctuation is less than 10 %. The spin lifetimes of the 
holes are distributed between 2.7 ms to 3.6 ms and the 
fluctuation of the lifetimes of holes is about 33 %. Anal- 
yses show that the spin lifetime fluctuation mainly comes 
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FIG. 3: (Color online) The electron- (a) and hole- (b) spin 
relaxation rates (l/Ti) of a pure InAs/GaAs QD as functions 
of temperature. QD base diameter=25 nm, height=3.5 nm. 
The black, red, and blue lines represent the spin relaxation 
rate due to the LADP, LAPZ, and TAPZ mechanisms, respec- 
tively. The green solid lines are the total relajcation rates. 



from fluctuations in the single-particle energy spacing 
(0.4 meV for electrons and 1.0 meV for holes) caused 
by the Ga distributions. 

Figure [3] depicts the spin-flip rates (green solid lines) 
for electrons and holes as functions temperature for a 
lens-shaped InAs/GaAs QD with diameter 13=25 nm, 
and height, /i=3.5 nm. The spin relaxation times for 
both electrons and holes increase quickly with decreas- 
ing temperature at low temperatures [9], because fewer 
phonons can take part in the scattering process. At lower 
temperatures, the holes may have extremely long life- 
times, for example, T=2 K, Ti=1.5 s, suggest that hole 
spins can be promising candidates for quantum comput- 
ing. The spin relaxation times T\ ~ 1/T^ for temper- 
atures >20 K, which is different from the one phonon 
process, where Ti ~ 1/T Q. We also compare the con- 
tributions of the three carrier-phonon interactions to the 
spin flip time. Among the three mechanisms, the contri- 
butions from the LAPZ (red dashed lines) and the TAPZ 
(blue dashes lines) are 2 to 3 orders of magnitudes weaker 
than that of the LADP (black dashed lines) from IK to 
room temperature (not shown). Below 1 K, the TAPZ 
has similar contribution to the spin flip rate to LADP. 

To summarize, we presented the flrst atomistic cal- 
culations of the multi-pho non-induced electron/hole spin 
relaxation in self- assembled Ini-^^Ga^^As/GaAs quantum 



dots at zero magnetic fleld. The calculated the spin life- 
times of the holes were in the range of a few millisec- 
onds to 20 milliseconds at 4.2 K, which is in excellent 
agreement with recent experimental results. We further 
studied the structural and alloy composition effects on 
the spin lifetimes, and found that the alloy QDs have 
much shorter hole-spin lifetimes. Compared with previ- 
ous methods, the current method has greater predictive 
power, and we expected it to play an important role in 
future studies on spin relaxations in quantum dots. 
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